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ABSTRACT 

Analytical expressions for the rates of formation and destruction of gravitationally 
bound systems are derived assuming that they are originated from primordial random- 
Gaussian density fluctuations. The resulting formulae reproduce the time derivative of 
the Press-Schechter mass function in a certain limit. Combining a theoretical model for 
the evolution of structures with the formation rate, we can make various cosmological 
predictions which are to be compared with observations. As an example to elucidate 
such applicability, we evaluate the contribution of clusters of galaxies to the cosmic X- 
ray background radiation. With the COBE normalization, we find that the significant 
fraction of the observed soft X-ray background is accounted for by clusters of galaxies 
in a cold dark matter universe with Jlo 0.2, Aq = 1 — r2o and h ^ 0.8. 

Key words: cosmology: theory - dark matter - clusters: formation - clusters: evolution 
- X-rays: sources. 



1 INTRODUCTION 

Formation, merging, and destruction history of dark mat- 
ter haloes of astrophysical objects is of fundamental impor- 
tance in describing the subsequent formation and evolution 
of luminous objects such as quasars, galaxies, and clusters 
of galaxies. In addition, early formation of structures may 
affect the growth of objects through the reionization of the 
entire universe (Sasaki, Takahara, & Suto 1993; Sasaki & 
Takahara 1994; Fukugita & Kawasaki 1994). 

Key quantities in this context are the rates of formation 
and destruction of gravitationally bound systems. They are 
defined as the comoving number density of bound systems 
of a given mass that are formed or destroyed in unit time 
at a given epoch. A number of authors have attempted to 
compute this quantity either numerically or analytically. A 
straightforward approach is to use a state-of-art numerical 
simulations (e.g., Evrard 1990; Cen & Ostriker 1992; Bryan 
et al. 1994; Suginohara 1994). Although the dynamical range 
and small-scale resolution of available codes are rapidly in- 
creasing, it is not clear at present to what extent they un- 
ambiguously describe small-scale dynamics. Besides it is not 
easy to grasp the underlying key physics from the results of 
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numerical simulations. It is therefore important to develop 
an analytical model on the basis of simple physics. 

The alternative approach in such a direction is closely 
related to the theory of Press & Schechter (1974; hereafter 
PS). PS derived a mass function, the number density as a 
function of mass and epoch, of virialized systems. The PS 
mass function, however, is not directly related to the rate 
of formation mentioned above. In fact, the time derivative 
of the PS mass function can become either positive or nega- 
tive. Thus it can not be simply interpreted as the rate of for- 
mation, but rather corresponds to the balance between the 
rates of formation and destruction (Cavaliere, Colafrancesco 
& Scaramella 1991; Blain & Longair 1993a, b). Although sev- 
eral authors have adopted the time derivative of the PS mass 
function to approximate the formation rate of bound sys- 
tems (e.g., Efstathiou & Rees 1988; Fukugita & Kawasaki 
1994), its validity should be examined carefully. Blain & 
Longair (1993a,b) and Sasaki (1994) correctly realized the 
above problem, and attempted to obtain the rates of forma- 
tion and destruction separately, assuming some empirical 
forms for the rates. In this paper, we derive a more general 
formula for each rate on the basis of merger probabilities dis- 
cussed by Bower (1991) and Lacey & Cole (1993; hereafter 
LC) (See also Kauffmann & White 1993). 

Our formalism is applicable to a number of problems re- 
garding the hierarchical clustering in the universe. In what 
follows we further apply the resultant formulae to the forma- 
tion of clusters of galaxies and calculate their contribution to 
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the observed cosmic X-ray background (XRB). The cluster 
contribution is in general supposed to be very small in the 
high energy band above a few keV, while it is likely to dom- 
inate in the soft energy band (Silk & Tarter 1973; Schaeffer 
& Silk 1988). It is, however, still uncertain to what extent 
the observed XRB can be contributed by emissions from 
clusters of galaxies. Several authors have employed the PS 
mass function to evaluate the cluster contribution (Evrard & 
Henry 1991; Blanchard et al. 1992; Burg, Cavaliere & Menci 
1993), and thus formation epochs and the subsequent evo- 
lution of virialized clusters are not well specified in their 
treatment. Here we present an approach which explicitly 
takes account of the formation epoch of clusters and their 
luminosity evolution. Our predictions are made specifically 
in cold dark matter (CDM) universes with/without the cos- 
mological constant whose fluctuation spectra are properly 
normalized by the COBE data (Sugiyama 1995). The re- 
sults are also compared with those of the hydrodynamical 
simulation and the ASCA observations. 

The plan of this paper is as follows. In Section 2, we out- 
line the PS formalism and discuss how the time derivative of 
its mass function is related to the formation and destruction 
rates of bound systems. Then using the merger probabilities 
of virialized haloes deduced by Bower (1991) and LC, we de- 
rive expressions for the rates of formation and destruction. 
These expressions are further employed to describe the dis- 
tribution of formation epoch of bound systems. In Section 
3, we apply our formulae and calculate the contribution of 
clusters of galaxies to the XRB. Finally Section 4 is devoted 
to our conclusions and further discussion. 



2 FORMATION AND DESTRUCTION RATES 
OF VIRIALIZED OBJECTS 

2.1 Press-Schechter mass function 

According to the hierarchical clustering model, small density 
fluctuations in the early universe grow via gravitational in- 
stability, become nonlinear, and eventually collapse to form 
bound virialized systems. For random Gaussian initial den- 
sity fluctuations, the comoving number density of spherical 
collapsed systems in the mass range M ~ M -I- dM at time 
t is evaluated as (Press & Schechter 1974; Bond et al. 1991) 
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where po is the mean comoving density of the universe, Sdt) 
is the critical density threshold for a spherical perturbation 
to collapse by the time t, and a{M) is the rms density fluc- 
tuation smoothed over a region of mass M. The latter two 
quantities are expressed respectively as 



<5c(t) 



and 



D{to) 

m 



Sc 

m 



{5c ~ 1.69 for ^0 = 1), (2) 



P{k,to)Wli{k)d''k, (3) 



a\M)^a\M,to) = j^ 

where D{t) is the linear growth factor normalized to unity 
at the present epoch to, fio is the cosmological density pa- 
rameter, P{k,to) is the present power spectrum of density 



fluctuations, and WM{k) is the Fourier transform of a real 
space window function whose volume contains mass M. 

Let us now consider the time derivative of the PS mass 
function (^): 
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The above expression is positive only for a{M) < 5c{t) since 
d5c{t)/dt < 0. This implies that Rps{M,t) is negative for 
low-mass objects for almost all realistic power spectra and 
window functions. Consider for instance scale-free spectra 
P{k) oc fc", which yield a'^{M) oc M-("+^'/^ where n > -3 
is physically plausible. The r.h.s. of equation (^ becomes 
negative for M < Mc{t) where Mc{t) is the critical mass at 
which cr(M) = Sc{t). This clearly indicates that RpsiM,t) 
should not be interpreted as the formation rate of bound 
objects. Rather it corresponds to the net rate of change in 
the number density. In other words, Rps {M, t) consists of the 
formation rate Rioim{M,t) at which objects of a given mass 
M are formed by mergers of smaller mass objects, and the 
destruction rate iidest {M, t) at which they are destroyed by 
coalescence to produce more massive systems. Thus equation 
should formally be written as 



(i) 



Rps{M,t) = i?form(A/,t) - i?dcst(M,t), 



(5) 



(Cavaliere et al. 1991; Blain & Longair 1993a, b). The neg- 
ative Rps{M,t) at M < Mc{t) is then ascribed to an ex- 
cess of the destruction rate over the formation rate in this 
mass range. On the other hand, the formation rate domi- 
nates at M > Mc{t) resulting in the positive Rps{M,t). In 
the case of scale-free spectra with n > —3, the critical mass 
Mc{t) increases monotonically with time because Mc{t) oc 
D(t)^^^"^^\ This is a generic feature in the 'bottom-up pic- 
ture' of hierarchical structure formation, in which small- 
scale structures form flrst and large-scale structures appear 
later. Our main task in this paper is to derive analytical ex- 
pressions for Rtc,Tin{M,t) and Rdcst {M,t). In so doing, we 
find that a more specific definition of formation and de- 
struction is essential in separating properly the two terms 
in equation (^. 



2.2 Conditional probabilities of merger process 

Merging is one of the most fundamental physical processes 
for the formation and destruction of virialized systems. LC 
developed a method to treat this process analytically follow- 
ing the idea of Bond et al. (1991). Their results are strictly 
justified only when one adopts the sharp fc-space filter- 
ing, i.e. the spherical top-hat filtering in fc-space, and when 
a{M) is a monotonically decreasing function of M. They 
found an analytical expression for the conditional probabil- 
ity that a point in a universe resides in an object of mass 
All ~ Ml + dMi at time ti provided it becomes part of a 
larger object of mass M2(> Mi) at later time t2{> ti): 



Pi{Mi,ti\M2,t2)dMi 
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where ai = a{Mi) and Sd = 5c{ti) for i = 1 and 2. The 
same expression was obtained in a more heuristic manner 
by Bower (1991). The reverse conditional probability that 
a point resides in an object of mass M2 ~ M2 + dM2 at 
time t2 provided it has been part of a smaller object of mass 
Mi(< M2) at ti(< t2) is 

P2{M2,t2\Mx,ti)dM2 
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Equations and (^) readily yield instantaneous transi- 
tion rates from a certain mass to another. Setting ti = t—At, 
t2 — t and A/2 = M in equation (P), and taking the limit 
At 0, one obtains 
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This quantity is interpreted as the rate at which an object 
of mass M{> Mi) is formed from an object of mass A/i ~ 
Ml +dMi in unit time at t. Similarly from equation one 
finds 
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Equation (^) corresponds to the rate at which an object of 
mass M{< M2) is incorporated into a larger object of mass 
A/'2 ~ A/2 + rfA/2 at time t. 

2.3 Formal procedure 

Now we can write down the rates of formation and destruc- 
tion of gravitationally bound systems using the transition 
rates in the last subsection. From equation (|8|), one formally 
defines the formation rate of an object of mass M at time t 
as 

X f'^ dPi(Mi ^ M-t) , 

/?form(A/,t) = / dMl '-LNps{M,t) 
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Equation (^) similarly defines the formal destruction rate: 
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Clearly these formal expressions diverge, reflecting the di- 
vergence of equations (b|) and (H) in the limit Mi M and 



A/2 M, respectively. Such a divergence would simply be 
related to the fact that the above formalism identifies the 
formation and destruction of objects even if the associate 
change of mass is infinitesimal. In reality, however, such a 
small change in mass hardly corresponds to the formation or 
destruction in a usual sense; rather it should be regarded as 
accretion or evolution of the same system. In order to incor- 
porate more realistic and physically reasonable description 
of these events, we should replace the above definitions of 
formation and destruction with more appropriate ones. Then 
the divergence associated in equations (^0^ and (|l^) would 
disappear naturally. 

One may argue that the divergence is ascribed to count- 
ing the transition of a certain object to practically the same 
one, by showing explicitly that equations (jl^) and (jl^) di- 
verge in a similar manner. In fact, one can prove that the 
divergences in the two rates exactly cancel out and their dif- 
ference reduces to the time derivative of the PS mass func- 
tion: 
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The above result also suggests that the divergence in equa- 
tions ( p^ ) and (jl^) can be removed by redefining the for- 
mation and destruction events properly. We will propose a 
realistic and reasonable prescription in the next subsection. 



2.4 Realistic prescription to compute Rfoim and 

/?dcst 

Since the divergence mentioned in the last subsection is 
originated from counting the transitions from an object to 
the essentially same one, a simplest remedy is to introduce 
thresholds in the mass of formation and destruction events. 
More specifically, we replace equations ( |lo| ) and (|ll| with 
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Physically the above procedure corresponds to assuming 
that an object of mass M keeps its identity as long as it 
stays between A/f and M^, and to defining the formation 
and destruction of the object only when the associate mass 
change is out of this range. Adopting such a prescription, 
the above formulae are integrated to give 
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Figure 1. Dependence of Ri^-^^^M , z; Mf) (solid curve), 
Rdast{M, z\ M(j) (daslied), and tlieir difference (dotted) on M/M{ 
or M^/M at 2 = in the standard CDM model (each quantity is 
normalized by \Rps{M , z)\); (a) M = IO^Mq, (b) M = lO^Afg, 
(c) M = 10^3 M0, (d) M = 10^5 M0. 
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where Z is defined as 
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ZiM,t;M,)^S.ity^ 2<.^(M).^(M.; 
and erfc(M) is the complimentary error function: 
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In the above prescription, it is important to choose 
physically reasonable values for Mf and M^- Admittedly it 
is inevitably artificial. We here propose to set Mf = M/2 
and Md = 2M as in LC mainly because of the following 
two reasons. One is that the identity of an object of mass 
M is clearly retained between M/2 and 2AI. If a new ob- 
ject of mass M were formed from a parent object of mass 
greater than M/2, the parent object would make up the 
dominant part of the new object. Therefore such a process 
should correspond to evolution of the parent object rather 
than formation of a new object. The same is true for the 
destruction of an object of mass M which is incorporated 
into an object of mass less than 2M. 

The other reason is more practical; equations ( [l5| ) and 
(^) are rather insensitive to the specific choice of Mf and 



Md at Mf « M/2 and Md ~ 2M, and their difference agrees 
fairly well with 7?ps(M, t). To see such features more clearly, 
we plot i?form(M, 2:; Mf ), _Rdest(M, 0; Md) and their differ- 
ence against M/Mf or Md/M in Fig. 1 (each quantity is 
normalized by |7?ps(M, z)\). Hereafter in this paper, a CDM 
power spectrum of Davis et al. (1985) is adopted and the 
mass variance is evaluated by the following approximation 
(White & Frenk 1991): 



1 - 0.4490(^0/1)° 



0.1 / ri 



/i-iMpc 



+ 0.6352(no/i) 



0.2 / n 



(19) 



where h = Ho/{100 km/sec/Mpc) is the dimensionless Hub- 
ble constant, and the top-hat filtering is implicitly assumed. 
The above approximation holds within 10% over the range 
0.013(^0/1)"^ < rf/(/i-iMpc) < 10(^0/1)"^ Equation ^ 
is normalized by the relation a{ri = 8/i~^Mpc) = , where 
b is the biasing parameter. In this section, we present results 
only in the standard CDM model (Qq = 1, /i = 0.5, and 
fe = 1) for an illustrative purpose, while we consider other 
CDM variants (e.g.. Table 1) in Section 3. 

Panels (a)-(d) in Fig. 1 are plotted for several masses at 
present redshift z = 0. The figure shows that i?form/|-Rps | 
and -Rdesst/I^psl are fairly insensitive to M/M; or M^/M for 
M/Mf, Md/M > 1.5, and that 7?form - i?dcst ~ -Rps within 
30% for M/M{,Mi/M < 3. Thus for practical purposes 
as well, the choice of Mf = M/2 and Md = 2M seems 
reasonable. 

Figure 2 illustrates 7?form(M, 2; M/2), i?dost(M, 2; 2M), 
their difference, and |i?ps(M, z)| versus M in the standard 
CDM model. Both panels exhibit that -Rform < -Rdost for 
small masses and iiform > ^dost for large masses. The critical 
mass Mc at which i?form = -Rdcst increases with time; Mc ~ 
1O^"M0 at 2 = 5, and Mc ^ IO^Mq at 2 = in this 
example. The difference between two rates is very close to 
RpQ except at M « Mc. The deviation simply arises from 
the fact that we set M{,Md 7^ M, but is not important in 
practice because its amplitude is several orders of magnitude 
smaller than those of i?form and 7?dcst • With these results in 
mind, we adopt Mf = M/2 and Md = 2M in the rest of 
this paper, and simply use i?form(M, t) and i?deat(M, t) to 
denote Riarm{M,t; M/2) and i?dost(M, t; 2M) respectively. 



2.5 Formation epoch distribution 

Before proceeding to the application of the results obtained 
in the previous subsections, let us discuss the formation 
epoch of virialized systems in some detail. First consider 
the probability that an object of mass M at time ti remains, 
without destructed, to have mass less than 2M at later time 
t2{> ti). We call this quantity the survival probability. It is 
evaluated by LC from equation (^) as 



Psurv(M,tl,t2) 



2M 
M 

2<5c2 



5ci 



P2{M2,t2\M,ti)dM2 
2(5c2(4l — 5c2) 



exp 



2(M) 
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Figure 2. Rates of formation i?forni(M, z, Aff) (solid curve) and destruction R^^^^(M, z, M^) (long-dashed) in the case of A/f = M/2 
and Afd = 2M, as a function of mass in the standard CDM model. Also plotted are the absolute values of their difference (short-dashed) 
and Rps{M, z) (dotted); (a) z = 0, (b) z = 5. 



<erfc[X(M,ti,t2)] - ierfc[F(M, ti, tz)], 



(20) 



where X{M,t\,t2) and Y{M,ti,t2) are defined respectively 



X{MMM) = 
and 

Y{M,ti,t2) = 



a^(2M)[gci - 25^2] + a^(M)^c2 
^2cr2(Af)a2(2A/)[a2(A'/) - a2(2M)] ' 



cr2(M)5c2-o-2(2A/)4i 



(21) 



(22) 



^2cr2(Af)a2(2Af)[a2(Af) - a2(2A/)] 

Using the above formulae and equation (^^, we write 
down the number density of bound systems which form with 
mass M ^ M -\- dM at time t{ ~ ti-\- dt{ and survive without 
destructed until later time t as 



F{M,ti,t)dMdtf 
1 



RiorUM,tf)dMdti X P,^,,{M,tt,t) 
'5c(tf) -2<5c(t) 



2tt y/a2{M/2)~a^{M) 
25,(t)[S4ti) - S4t)] 



2 - 



X exp 



a2(Af) 
erfc[y(Af, U,t)] 



dSc , , 



Sc{t) 
erk[X{M, tf,t)] 

Nps{M, tt)dMdtf. (23) 



Given a mass M and time t, equation ( |23[ ) is regarded as the 
distribution function of formation epoch t{. In Fig. 3 we plot 
F{M, Z{, z = 0) versus Zf for several masses in the standard 
CDM model. A common feature is that f as a function 
of zt remains nearly constant up to a certain redshift and 



then starts to decline. The epoch at which F drops by an 
order of magnitude is 2f « 6 for M = 10^ A/©, Zf ~ 4 for 
M = 10"A^o, Zf X 2 toT M = IO^^Mq, and zt ^ 0.6 for 
M = 10^^ Mq. These are interpreted as typical epochs for 
the formation of bound systems of respective masses which 
exist at z — 0. 



3 X-RAY BACKGROUND RADIATION 
CONTRIBUTED FROM CLUSTERS OF 
GALAXIES 

As a simple and important application of our formulae, we 
calculate X-ray emission from clusters of galaxies and their 
contribution to the observed XRB. The formation rate of 
virialized systems plays an essential role in specifying the 
epoch of formation of clusters and in explicitly taking ac- 
count of their luminosity evolution. 



3.1 Method of calculation 

3.1.1 Cluster contribution to the XRB spectra 

Clusters of galaxies are filled with diffuse hot plasma which 
produces copious X-ray emission. If the total power liberated 
per unit comoving volume and unit frequency at redshift 
z is denoted by J,^{z), the background X-ray intensity 
observed at present at frequency uq is described as 
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M = 1.0E+09 [Mg,,}^ 
M = 1.0E+11 [M3JE 

M = 1.0E+13[M3,,n 
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Figure 3. Distribution of formation epoch F{M,Zf,z) in the 
standard CDM model. The differential number densities of bound 
objects that are formed at redshift Z{ and survive until z = are 
plotted against Zf for masses M = 10^ Mq (solid curve), 10^^ Mq 
(dotted), lO^^M© (short-dashed) and lO^^M© (long-dashed). 



c f°° 
4^70 (1 



J, 



i^o(l-l-z) 



(z) dz 



Jo (1 + z)y'{floz + 1)(1 + zf - Ao2(2 + z) 



.(24) 



where c is the speed of light, and Ao is the dimensionless 
cosmological constant. In our model, the luminosity at fre- 
quency V oi a. cluster, L^(M, 2f,z), is assumed to be deter- 
mined by the mass M, the epoch of its formation Zf, and the 
epoch z at which the observed X-ray is emitted (see eq.|34| 
below). Then Jv{z) is given by 



dM 



dzcF(M, Zf, z)L^{M, z{, z) 



(25) 



where F{M,Zf,z) is the formation epoch distribution 
(eq. H). 

In contrast to equation (Pq), previous authors calcu- 
lated Jv{z) using the mass function of clusters directly as 
follows (Evrard & Henry 1991; Blanchard et al. 1992; Burg 
et al 1993): 

poo 

J4z)= dMN(M,z)L,(M,zi = z,z), (26) 
Jo 

where the mass function N{M, z) is often taken as the PS 
formula (eq. [Q|). Since the mass function does not specify 
the epoch of formation of virialized clusters, zi is usually 
set to be equal to the epoch of emission z. Such a treat- 
ment, however, is strictly justified only if Zi is very close to 



z, and will become inappropriate otherwise. The properties 
of a virialized system such as radius and temperature are 
determined by Zf not z (see Section j.l.S ). Thus the lumi- 
nosity evolution described in equation (p6D is not realistic. In 
fact one frequently encounters similar situations where the 
formation epoch is essential in incorporating the evolution 
of individual astronomical objects and making cosmologi- 
cal predictions. This is one of the strong advantages of our 
current refinement of the original PS theory. 



3.1.2 X-ray emission from clusters 

The total gravitating mass of clusters is dominated by dis- 
sipationless dark matter, which is likely to be relaxed im- 
mediately after the collapse of halo and attain the virial 
equilibrium. In a theoretical treatment for the collapse of 
spherical perturbations, the radius of a virialized cluster in 
the Einstein-de Sitter (Qq ~ 1) universe is estimated as a 
half of maximum turn-around radius and given by 



(GM) 



1/3 



(37rffo)2/3 1 + Z{ 
1.69(1 + 2f)" 



M 



IQi^Mo 



1/3 



/i"^/^ Mpc, (27) 



where M is the total mass of the cluster, and Zf is the redshift 
of its formation. The virial temperature of the system Tvir 
is then calculated as 
_ GMfj,mp 



3fcB»'v 



M 



\^ 1015 Me 



2/3 



l2/3 



K, (28) 



where G is the gravitational constant, ks is the Boltzmann 
constant, nip is the proton mass, and /i is the mean molec- 
ular weight. The cases of fio < 1 universes are discussed 
in Appendix A. Hereafter we assume that the intracluster 
gas is fully ionized with the primordial abundance Y — 0.24, 
where Y is the helium fraction by weight, and thus /i — 0.59. 
Note that once a model for the background universe is spec- 
ified, properties of virialized clusters are determined by two 
parameters; mass M and the formation epoch z{ which is 
not available in the conventional PS theory. 

The gas component of clusters generates the X-ray 
continuum emission mainly due to thermal bremsstrahlung 
whose emissivity in CGS units is (e.g., Rybicki & Lightman 
1979) 



= 6.84 X 10" 



-1/2 



exp 



Y 
hi' 



9{Tg: 



BJgas 



ergs-icm-=*Hz-f29) 



where h is the Planck constant, Tgas is the gas temperature, 
and Tic is the electron density. We adopted an approximation 
g{T,iy) = 0.9(/w'/feBr)"°-^ for the velocity -averaged Gaunt 
factor following Evrard & Henry (1991). 

The X-ray emission described by equation ( p9| ) is es- 
sentially determined by the temperature and density pro- 
files of intracluster gas. Simple models for the intraclus- 
ter gas often assume that the gas is isothermal and its 
distribution is spherically symmetric. The observed X-ray 
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spectra and galaxy velocity dispersions further suggest that 
Tvir/Tgas ^ 1 (Mushotzky 1984). Thus we set 



Lboi(M, zc, z) oc M^/3(l + ziY^^il + z) 



Tgas = Tvir(Af, 2f). 



(30) 



For the density profile, the observed X-ray image of surface 
brightness is found to agree well with (Jones & Forman 1984) 

-1 



Pgas(r) = Pg 



l+(^ 



(31) 



where fP^^^ is the central gas density, r is the distance from 
the cluster centre, and Tc is the core radius. The above mod- 
eling is slightly different from the conventional isothermal /3- 
model, in which the temperature and density profiles of intr- 
acluster gas are described by a single parameter /3. Although 
our parametrization seems inconsistent with the /3 model, 
it can be reconciled if the galaxy number density scales as 
[1 + (rjr^fY^ unlike the King model [1 -h {rjr^fY^I'^. In 
fact, this is qualitatively indicated from some observations 
(BahcaU & Lubin 1994). 

As for the core radius rc, we still lack understandings 
of its physical meaning, and are unable to predict how it 
depends on M and Zi. Thus we adopt a simple self-similar 
model in which rc is proportional to rvir: 



0.15/i~^Mpc X 
0.15(1 



rvir(Af, Zi) 



Zi) 



rvir (1015 Mo, = 0) 

1/3 

/i"^Mpc, 



M 
1015 M( 



(32) 



where the normalization is chosen to match the Ginga obser- 
vations (Hatsukade 1989), and the second equality assumes 
the Einstein-de Sitter universe. Then the central gas density 

Pi 



is fixed by 



(33) 



where Q.b is the baryonic density parameter. 

In order to describe realistic X-ray emission from clus- 
ters, we also need to take account of intrinsic luminosity 
evolution. As the effects of merger are already factored out 
in equation (psj), quiet gas accretion and the subsequent 
changes in temperature and density profiles are here to be 
incorporated. These processes, however, are still highly un- 
certain. So we attempt to model their effects into a factor 
^{zi, z) and express the luminosity as 



L„{M, zi,z) = 7(zf, z) 



el/ {M, Zi; r) Anr^dr, 



(34) 



where L^{M, zt, z) is the luminosity that a cluster formed 
with mass M at redshift Zi would have at a later epoch z. 
Given a large degree of uncertainty, 7(2:5, z) is assumed here 
to take a power law form: 



(35) 



Current numerical simulations seem to indicate p ^ 1 (e.g., 
Evrard 1990; Suginohara 1995; Navarro, Frenk & White 
1995), and we examine two limiting cases of no intrinsic 
evolution, p = 0, and of very strong evolution, p — 1. 

Integrating equation (134) over whole frequency, one 
finds a simple scaling law for the bolometric luminosity of 
the form: 



(36) 



where the parameter q is constant and equal to 7/2 in the 
Einstein-de Sitter universe, while it varies with Zi and less 
than 7/2 in low flo universes; e.g. 1.8 ^ q ^ 3.5 for fio ~ 0.2 
and Ao = 0.8. 

3.1.3 X-ray temperature and luminosity functions 



In Section 3.1.1, we have emphasized the difference between 



our formulation and the conventional PS approach. For more 
direct comparison of these two procedures, we also calculate 
the X-ray temperature and luminosity functions of clusters 
of galaxies. 

In our model, the X-ray temperature function at red- 
shift z is expressed using equation (ESI) as 



NT{T,z)dT^ 



' , ,dM 
dziF{M,Zi,z) — 



dT, (37) 



M=M{T,zi) 



where M{T, Zi) stands for the mass which gives the gas tem- 
perature T if formed at redshift Zf (e.g., eq. |^). Here we 
have neglected the intrinsic evolution of gas temperature. 
On the other hand, from the PS mass function, one obtains 



NT{T,z)dT = Nps{M,z) ^ 



dT. (38) 

M=M{T,Zf=z) 

Note that Zi is set to be equal to 2 in the above equation. 

Similarly, the X-ray luminosity function at redshift 2 in 
each approach is 



N^{L,z)dL ■ 



and 



dziF(M, zt,z) — 



dL, (39) 



M=M{L,Z{,z) 



Ni,{L,z)dL = Nps{M,z) 



dM 



dL, 



M=M(L,zf-z,z) 



(40) 



where M{L, Zi, 2) is solved from the integral of equation ( |34[ ) 
over frequency. 

3.2 Comparison with hydrodynamical simulation 

It is meaningful to check the validity and limitations of our 
modeling through the comparison with other classes of ap- 
proach. Fortunately a series of three-dimensional hydrody- 
namic simulations in a CDM model universe have been de- 
veloped to study the evolution of X-ray clusters (Cen et al. 
1990; Cen & Ostriker 1992, 1993, 1994; Kang et al. 1994; 
Bryan et al. 1994). In this subsection, therefore, we com- 
pare the results of our analytical approach with those of 
numerical simulations. 

Figure 4 illustrates the temperature functions evalu- 
ated analytically by equations (^) and ( p^ ) respectively, 
in comparison with the results of the simulation by Kang 
et al. (1994) and the observations (Henry & Arnaud 1991). 
The cosmological parameters are the same as in Kang et 
al. (1994); D,o = 1, h = 0.5, D.b = 0.06 and b = 0.95. The 
indicated error bars of the simulation data correspond to sta- 
tistical ones only. The simulation data are kindly provided 
by Renyue Cen, and note that the original plots in Figures 
1 to 5 of Kang et al. (1994) should be shifted down by a 
factor of In 10 due to an error in their plotting procedure 
(Cen, private communication). 



8 T. Kitayama and Y. Suto 




10"' 1 10' 

T [keV] 



Figure 4. X-ray cluster temperature functions at redshift 2 = 
in the CDM model; predictions of our model (eq. [p7|) (solid 
curve), of the PS approach (eq. |pq]) (dotted curve), and of the 
simulation by Kang et al. (1994) (crosses). The cosmological pa- 
rameters are fixed in each as SIq = 1, h = 0.5, Ob = 0.06 and 
b = 0.95. The error bars of the simulation data only indicate 
the statistical error. The squared region shows the error box of 
observations (Henry & Arnaud 1991). 




10"^ 10"^ 10"' 1 io' 
L,„,[10''h"' erg/s] 

Figure 5. X-ray cluster bolometric luminosity functions at red- 
shift 2 = in the CDM model; predictions of our model (eq. B? ) 
with the parameter p = (solid curve)and p = 1 (dotted curve), 
of the PS formalism (eq. [M) (dashed curve), and of the simula- 
tion by Kang et al. (19941 (crosses). The cosmological parameters 
are the same as in Fig. W The error bars of the simulation data 
are only statistical. The squared region shows the error box of 
observations (Henry & Arnaud 1991). 



The analytical results in Fig. 4 show a reasonable agree- 
ment with that of the numerical simulation above T ~ IkeV. 
In fact, the simulation data are subject to much larger er- 
ror below this temperature due to its spatial resolution of 
~ 0.75/i"^Mpc (Kang et al. 1994). They are also likely to be 
afTected by the limited sample volume size at T > lOkeV. 
The difference between our method and the previous PS 
approach is in general rather insignificant in this figure, 
because the temperature of a virialized system depends 
only weakly on Zf as oc (1 + Zf) in the fio = 1 universe 
(eq. [^). Such a feature is especially apparent at high tem- 
peratures, where most clusters that contribute are large and 
formed only at low redshifts. The difference becomes rela- 
tively larger at low temperatures where small clusters make 
dominant contributions. 

The results for the bolometric luminosity functions are 
plotted in Fig. 5. Neither the simulation nor our analyti- 
cal approach takes account of physical cooling or heating. 
This is justified in discussing the overall properties of clus- 
ters of galaxies except at the core. In contrast to Fig. 4, 
the difference between our method and the conventional PS 
approach is more significant, because the bolometric lumi- 
nosity depends much more strongly on Zi as oc (1 + ZfY^"^ 



(eq. Igg). The deviations from the simulation data are also 
larger and our changes to the PS formalism does not rec- 
oncile the discrepancy. However, as remarked by Cen (1992) 
and Kang et al. (1994), the numerical procedures tend to sys- 
tematically underestimate the predicted luminosity, perhaps 
by a factor of ~ 2. The underestimation would be severer 
at low luminosities where the effects of the spatial resolu- 
tion becomes more important. In addition, the simulation 
clearly lacks the sample volume size for the statistically re- 
liable predictions at Lx ^ 10'''*/i~^ erg/s where most of the 
observed data are available (Henry & Arnaud 1991). Thus 
the analytical results are not only qualitatively consistent 
with the numerical simulation, but also applicable in much 
wider dynamical range. 

Figure 6 shows the predictions of the XRB spectra, to- 
gether with the ASCA observations. The analytical results 
agree well with the simulation around E ^ IQ ke'V but pre- 
dict larger intensity at i5 < 2 ke'V. This is reasonable if one 
considers the numerical effects including the spatial reso- 
lution of the simulation. As expected from the features in 
Fig. 5, our model systematically predicts higher intensity 
than the original PS formalism, and the intensity increases 
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Figure 6. The XRB spectra contributed from clusters of galaxies; 
predictions of our model (eq. [p5|) with the parameter p = (solid 
curve) and p = 1 (dotted curve), of the PS formalism (eq. p6[ ) 
(dashed curve), and of the simulation by Kang et al. (1994) (tri- 
angles). The cosmological parameters are the same as in Fig. W. 
Also plotted are the observed data from ASCA with statistical 
errors (Gendreau et al. 1994). 



as the intrinsic evolution indicated by the parameter p be- 
comes stronger. 

As a matter of fact, our method and numerical simu- 
lations are complementary to each other in predicting the 
X-ray properties of clusters and their contribution to the 
XRB spectra. The analytical approach basically has two ad- 
vantages; one is that it is much easier to make systematic 
parameter surveys. Therefore it would give a better insight 
into how the predictions depend on the distribution of intr- 
acluster gas, the luminosity evolution, and the geometry of 
the background universe. The other is that the results are 
largely free from the small-scale resolution or the limited 
volume size that may affect the numerical results. Alterna- 
tively, our present predictions do not properly take account 
of the thermal evolution of clusters induced by the interac- 
tion between gas and radiation, which will be incorporated 
in due course. In addition, we have assumed a rather sim- 
plified gas density profile as well as a spherical distribution 
of dark matter. 



Table 1. Models for the background universe 



Model 




Ao 


h 




b 


E5 


1 





0.5 


0.05 


0.71 


E8 


1 





0.8 


0.02 


0.43 


05 


0.2 





0.5 


0.05 


6.7 


08 


0.2 





0.8 


0.02 


2.8 


L5 


0.2 


0.8 


0.5 


0.05 


2.2 


L8 


0.2 


0.8 


0.8 


0.02 


1.0 



3.3 Predictions for various cosmological models 

In this subsection, we present the XRB spectra to constrain 
the viable sets of cosmological parameters in a CDM uni- 
verse. Specifically we consider six CDM models listed in 
Table 1. Values of the Hubble parameter h are taken as 
either 0.5 or 0.8, where the latter is strongly motivated 
from the observation of Cepheids in the Virgo cluster by 
the Hubble Space Telescope (Freedman et al. 1994). In each 
model, the baryon density of the universe is assumed to be 
Qsh'^ = 0.0125, which is consistent with primordial nucle- 
osynthesis (e.g., Walker et al. 1991). The biasing parame- 
ter b is given by normalizing the amplitude of the fluctua- 
tion spectrum by the COBEDMK two year data (Sugiyama 
1995). The critical overdensity for the spherical collapse Sc 
is taken as Sc = 1.69 in models E5, E8, L5 and L8, and 
5c = 1.64 in 05 and 08 (Lilje 1992; LC). 

Figures 7, 8 and 9 illustrate the calculated XRB spec- 
tra together with the ASCA observations (Gendreau et al. 
1994). There are a few points that need to be kept in mind. 
Firstly, line emissions or absorptions from metals are omit- 
ted in our calculation. This may lead to underestimation of 
XRB intensity below 1 keV (Cen et al. 1995). Secondly, a 
steep rise in the observed intensity in the soft X-ray band 
may be partly due to Galactic emission (McCammon & 
Sanders 1990; Snowden et al. 1990). Finally, it has been 
suggested from the high-resolution image of RO SAT that at 
least 30% of the XRB at 1 keV is contributed from quasars 
(Shanks et al. 1991). 

Having considered the above points, we see that models 
E5 and E8 overproduce the soft X-ray intensity (Fig. 7). 
While this result was pointed out earlier by Burg et al. 
(1993) on the basis of the conventional PS approach, such 
a feature is more apparent in our model which takes ac- 
count of the formation rate. The overproduction could be re- 
moved if the clusters underwent considerable anti-evolution 
against the growth of the gravitational potential (Burg et 
al. 1993). Alternative ways to avoid such discrepancy are to 
consider a small value for the Hubble constant; e.g., h < 0.3, 
or to abandon the COBE normalization and introduce much 
smaller fluctuation amplitude, both of which are highly arti- 
ficial even if not completely excluded. Therefore the Slo = 1 
CDM universes have a serious difficulty of overproducing the 
XRB. 

On the other hand. Fig. 8 shows that the predicted in- 
tensity in models 05 and 08 is significantly smaller and that 
clusters of galaxies do not basically contribute to the XRB 
in these models. As a matter of fact, this is mainly because 
of large values of the biasing parameter b imposed from the 
COBE data. Clusters of galaxies are thus rejected as a can- 
didate for the origin of the XRB in these open models. 

Noticeably, Fig. 9 suggests that the observed XRB at 
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Figure 7. The XRB spectra contributed by clusters of galaxies in the Einstein-de Sitter universe; (a) E5, (b) E8. Lines show predictions 
of our model with the parameter p = (solid) and p = 1 (dotted), and of the original PS formalism (dashed). Also plotted are the ASCA 
observations from Gendreau et al. (1994). 
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Figure 8. Same as in Fig. 7 except that the background universe is open (Qq = 0.2, Ao =0); (a) 05, (b) 08. 
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Figure 9. Same as in Fig. 7 except that the background universe is spatially flat with the cosmological constant (f2o = 0.2, Xq = 0.8); 
(a) L5, (b) L8. 



E < O.SkeV attains a dominant contribution from clusters 
of galaxies in model L8. We should, however, be careful be- 
cause the observation in this band is most sensitive to the 
contamination due to Galactic absorptions and emissions. 
Nevertheless it is interesting to see a possibility that a part 
of the XRB can be naturally explained by clusters of galaxies 
in one of promising cosmological models. 



4 CONCLUSIONS 

We have extended the theory of Press & Schechter (1974) 
and the formulation of Bower (1991) and Lacey & Cole 
(1993) to derive expressions for the rates of formation and 
destruction of gravitationally bound systems. Although the 
formal expressions of these quantities diverge, they properly 
reproduce the time derivative of the PS mass function. The 
divergence is not unphysieal but can be interpreted as being 
originated from the somewhat improper definition of forma- 
tion and destruction of bound objects. We have proposed a 
phenomenological prescription to remove such divergences 
by specifying a realistic mass range in which the identity of 
an object is to be maintained. The resulting expressions are 
practically insensitive to the choice of the threshold masses. 
On the basis of this prescription, we have obtained a distri- 
bution function of formation epochs of bound systems. This 
quantity is of great cosmological importance because it pro- 
vides a useful theoretical tool in describing the subsequent 
evolution of the systems. 

In order to exhibit the applicability of the above formal- 
ism, we have evaluated the contribution of clusters of galax- 



ies to the XRB, as well as the X-ray cluster temperature 
and luminosity functions. Changes from the previous PS for- 
malism are obvious especially when the cluster evolution is 
strong and the X-ray properties depend largely on its forma- 
tion epoch. The results of our analytical procedure, although 
fairly phenomenological in modelling clusters' X-ray gas, are 
reasonably supported by those of the three-dimensional hy- 
drodynamical simulation by Kang et al. (1994). Thus these 
two approaches are complementary in investigating the evo- 
lution history of virialized structures. 

We have shown that the significant fraction of the ob- 
served soft XRB can be accounted for by clusters of galax- 
ies in a spatially fiat CDM universe with ilo ~ 0.2 and 
Ao = 1 — fio with the COBE DMR two year normalization. 
The f2o = 1 CDM model tends to overproduce the observed 
X-ray intensity. On the other hand, the cluster contribution 
becomes negligibly small in an open CDM universe with 
0.0 ~ 0.2, in which case observed XRB should be explained 
entirely by other classes of sources such as quasars. Although 
our overall conclusion in the f2o = 1 case agrees with Burg 
et al. (1993), we derived new constraints on the other mod- 
els from the XRB in a consistent and clear manner. Note 
that the model with f2o ~ 0.2, Ao = 1 — Oq and h ^ 0.8 is a 
specific example of the most successful cosmological models 
explaining the cosmic age, two-point correlation functions 
of galaxies and clusters, and small-scale velocity dispersions 
(Efstathiou et al. 1990; Suginohara & Suto 1991; Suto 1993; 
Watanabe, Matsubara & Suto 1994; Bahcall 1994). In this 
context it is interesting to note that the same model can 
simultaneously account for a major fraction of the observed 
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XRB in the soft energy band. Furthermore, the microwave 
background distortion and the resulting anisotropy due to 
the distant clusters in CDM models are below the current 
limit (Makino & Suto 1993). 

Obviously there remain several areas for further theo- 
retical work. With regard to the XRB predictions, the great- 
est uncertainty probably lies in the evolution of intracluster 
medium. Although we have assumed that it basically traces 
the growth of the gravitational potential, this assumption 
is still an open question and different evolutionary models 
would alter the predicted luminosity functions and the XRB 
spectra. The contamination of the observed spectra by the 
Galactic emission and absorption is also uncertain. More- 
over, we have omitted the line emissions which may affect 
the soft band intensity. Therefore we did not aim to rigor- 
ously reproduce the XRB spectra from theory, but rather 
to examine the possibility that clusters of galaxies play a 
crucial role in the origin of the XRB. 

As to the formalism, on the other hand, our method 
contains a formal divergence for which somewhat artificial 
prescription has been proposed. In addition, our argument 
itself has further room for improvement; for instance, if Mi 
in equation ( |l^ ) does not correspond to the largest parent 
of the final object M, the rate of formation thereby defined 
possibly overcounts the true rate. This problem is in fact 
ascribed to a fundamental limitation of the current formal- 
ism which does not take into account the history of every 
single object involved in a merger process. Clearly the valid- 
ity of our treatment needs to be examined through a direct 
comparison with other classes of approach such as N-body 
simulations. Other ways of formulation are also desirable. 
In this context we note that LC proposed a differential dis- 
tribution function of halo formation times based on what 
they called 'the halo counting argument' (see their §2.5.2). 
As notified by LC, however, their definition leads to a slight 
self-inconsistency of predicting a negative probability den- 
sity depending on the slope of the fiuctuation spectrum. As 
a matter of fact, we compared their result with ours and ba- 
sically confirmed that both agree well in the mass range of 
astrophysical interest. For practical purposes, therefore, our 
present formulae are expected to be applicable to a num- 
ber of further problems in cosmology such as cluster lumi- 
nosity functions in various cosmological models, formation 
of high redshift objects and reionization history of the uni- 
verse. These issues, together with a quantitative comparison 
with LC, will be described elsewhere (Kitayama & Suto, in 
preparation). 



ACKNOWLEDGEMENTS 

We thank Shin Sasaki, Katsuhiko Sato, Tatsushi Sugino- 
hara, Hajime Susa, and Takahiro Tanaka for helpful dis- 
cussions and comments. We are grateful to Renyue Cen, 
Keith Gendreau and Naoshi Sugiyama for providing the rel- 
evant results of their current work which are adopted in 
the present paper. TK acknowledges the fellowship from the 
Nippon Telephone and Telegram Co. . This research is sup- 
ported in part by the Grants-in-Aid by the Ministry of Edu- 
cation, Science and Culture of Japan (05640312, 06233209, 
07740183, 07CE2002). 



REFERENCES 

Bahcall N. A., 1994, preprint POP-606 
Bahcall N. A., Lubin L. M., 1994, ApJ, 426, 513 
Blain A. W., Longair M. S., 1993a, MNRAS, 264, 509 
Blain A. W., Longair M. S., 1993b, MNRAS, 265, L21 
Blanchard A., Wachter K., Evrard A. E., Silk J., 1992, ApJ, 391, 
1 

Bond J. R., Cole S., Efstathiou G., Kaiser N., 1991, ApJ, 379, 
440 

Bower R. J., 1991, MNRAS, 248, 332 

Bryan G. L., Cen R., Norman M. L., Ostriker J. P., Stone J. M., 

1994, ApJ, 428, 405 
Burg R., Cavaliere A., Menci N., 1993, ApJ, 404, L55 
Cavaliere A., Colafrancesco S., Scaramella R., 1991, ApJ, 380, 15 
Cen R., 1992, ApJS, 78, 341 

Cen R., Jameson A., Liu F., Ostriker J. P., 1990, ApJ, 362, L41 
Cen R., Kang H., Ostriker J. P., Ryu D., 1995, submitted to ApJ 
Cen R., Ostriker J. P., 1992, ApJ, 393, 22 
Cen R., Ostriker J. P., 1993, ApJ, 417, 404 
Cen R., Ostriker J. P., 1994, ApJ, 429, 4 

Davis M., Efstathiou G., Frenk C, White S., 1985, ApJ, 292, 371 
Efstathiou G., Rees M. J., 1988, MNRAS, 230, 5p 
Efstathiou G., Sutherland W. J., Maddox S. J., 1990, Nature, 
348, 705 

Evrard A. E., 1990, ApJ, 363, 349 
Evrard A. E., Henry J. P., 1991, ApJ, 383, 95 
Freedman et al. 1994, Nature, 371, 757 
Fukugita M., Kawasaki M., 1994, MNRAS, 269, 563 
Gendreau K. C, et al., 1994, in Makino F., Ohashi T., eds. New 
Horizon of X-ray Astronomy —First Results from ASCA. Uni- 
versal Academy Press, Tokyo, p. 365 
Hatsukade I., 1989, Osaka University PhD thesis (unpublished) 
Henry J. P., Arnaud K. A., 1991, ApJ, 372, 410 
Jones C, Forman W., 1984, ApJ, 276, 38 
Kang H., Cen R., Ostriker J. P., Ryu D., 1994, ApJ, 428, 1 
Kauffmann G., White S. D. M., 1993, MNRAS, 261, 921 
Lacey C. G., Cole S., 1993, MNRAS, 262, 627 (LC) 
Lahav O, Lilje P. B., Primack J. R., Rees M. J., 1991, MNRAS, 

251, 128 
Lilje P. B., 1992, ApJ, 386, L33 
Makino N., Suto Y., 1993, ApJ, 405, 1 

McCammon D., Sanders W. T., 1990, Annu. Rev. Astron. Astro- 

phys., 1990, 28, 657 
Mushotzky R. F., 1984, PhyScr, T7, 157 

Navarro, J.F., Frenk, C.S., White, S.D.M., 1995, MNRAS, 275, 
720 

Peebles P. J. E., 1980, The Large-Scale Structure of the Universe, 
Princeton Univ. Press, Princeton 

Press W. H., Schechter P., 1974, ApJ, 187, 425 (PS) 

Rybicki G. B., Lightman A. P., 1979, Radiative Processes in As- 
trophysics. Wiley, New York 

Sarazin C. L., 1988, X-ray Emission from Clusters of Galaxies. 
Cambridge Univ. Press, Cambridge 

Sasaki S., 1994, PAS J, 46, 427 

Sasaki S., Takahara F., 1994, 91, 699 

Sasaki S., Takahara F., Suto, Y., 1993, Prog. Theor. Phys., 90, 
85 

Schaeffer R., Silk J., 1988, ApJ, 332, 1 

Shanks T., Georgantopoulos I., Stewart G. C, Pounds K. A., 

Boyle B. J., Griffiths R. E., 1991, Nature, 353, 315 
Silk J., Tarter J., 1973, ApJ, 183, 387 

Snowden S. L., Cox D. P., McCammon D., Sanders W. T., 1990, 

ApJ, 354, 211 
Suginohara T., 1994, PASJ, 46, 441 
Suginohara T., 1995, personal communication 
Suginohara T., Suto Y., 1991, PASJ, 43, L17 
Sugiyama N., 1995, submitted to ApJS 



Formation rate of gravitational structures 



Suto Y., 1993, Prog.Theor.Phys., 90, 1173 

Walker T. P., et al., 1991, ApJ, 376, 51 

Watanabe T., Matsubara T., Suto Y. 1994, ApJ, 432, 17 

White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 



APPENDIX A: RADIUS OF VIRIALIZED 
CLUSTERS IN LOW Qo UNIVERSES 

In an open universe (Qq < 1, Ao =0), the effective radius 
that a spherical perturbation attains when it virializes is 
estimated as half of its maximum radius at turn-around. If 
a virialized system is formed with mass M at redshift Zf, 
one can readily show that the radius r^ii is given by 



I- no { 47r//o ' ^^^^ 

where rjf is expressed in terms of the epoch of collapse Zf as 

2(1-^0) 



Tji = arccosh 



1 



(A2) 



In the presence of the cosmological constant Ao, the 
dynamical motion of a spherical perturbation is modified 
and the radius rvir of a virialized cluster is no longer half 
of the turn-around radius rta.. In fact, the ratio of these 
two quantities is shown to obey the following cubic equation 
(Lahav et al. 1991): 

2x(— )'-(2 + X) ( — )+! = 0, (A3) 
where 

X^'^rL (A4) 

The condition for a shell to turn around is x < 1. Lahav et al. 
(1991) found a following approximated solution to equation 
®: 

rvir ^ 2 - X 
rta 4 - X ' 

For Ao ~ 0, the above result reproduces rvir/rta = 1/2. 

Now the maximum turn-around radius rta is no longer 
expressed as a simple analytical form if Ao is nonzero. For 
our present purpose, however, we simply assume that the 
presence of Ao affects the dynamics of the perturbation only 
through equation ( |A3| ) and it does not alter the value of rta 
itself. Then in a spatially flat universe (f2o <1,A = 1 — f^o) 
which we discussed in the present paper, rta is written down 
as 

--^(^^^^^^( 3.^0^ )^'^' 
where Q is related to Zf via 



(A5) 



({ = arcsinh 



(l + zf)3 



(A7) 



